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1 Introduction 


When it comes to the question of the formation of a new stable phase from a 
supersaturated one, the frequency with which nuclei form is of primary interest. 
This view underlies almost all work that has dealt with the kinetics of the formation 
of a new phase, partly theoretically and partly experimentally. A nucleus was 
understood to be a small, spontaneously formed preferred droplet or crystal from 
which the visible particles of the new phase form as they grow. 

Until recently, there were no quantitative approaches to the size of the nucleus and 
its possible connection to the supersaturation pressure. The nucleus received a 
physical meaning through VOLMER and WEBER’s definition!, which understands by 
nucleus a droplet or crystal that has the same vapor pressure as the supersaturation 
pressure. This nucleus is in equilibrium with the supersaturated phase. This also 
determines the size of the nucleus, and a very specific amount of work ¢€ is required 
to produce this characteristic droplet inside the homogeneous vapor. This work 
was already calculated by GIBBS using thermodynamic methods and, in addition to 
the supersaturation pressure, depends on the capillary forces between the nucleus 
and the vapor. It is 


1 
e=-F-o, 
go 


where -¥ is the surface of the nucleus and a its surface tension. 

VOLMER and WEBER were the first to make an approach to the number of 
nuclei present in a supersaturated vapor in which, based on EINSTEIN’s general 
thermodynamic considerations’, they estimated the stationary number of nuclei on 


'VoLMER and WEBER, Keimbildung in iibersattigten Gebilden. Zeitschr. f. physik. Chemie 119, 
277 (1926). 
2A. EINSTEIN, Ann. d. Physik 33, 1275 (1910) 
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the order of 
n=C-e ur. 


Furthermore, in the cited work by VOLMER and WEBER it is assumed that the 
number of nuclei that grow into visible structures per second is proportional to the 
stationary number of nuclei. The model is the ARRHENIUS theory of the reaction 
rate, where the activation heat plays the same role as the characteristic work here: 
é. 

The aim here is to show that the above approach can be precisely justified on the 
basis of simple kinetic ideas and that the number of droplets that grow into a 
visible structure can be calculated from the stationary number of nuclei.?. 


2 Theoretical Considerations 


We want to limit ourselves to the case of supersaturated vapor, since the conditions 
are simplest here. 

Our basic assumption is that as a result of the collisions of the molecules in the 
vapor space, some molecules can combine, and that the particle number n or the 
radius x of the resulting droplet is subject to irregular size changes as a result of 
natural fluctuations. The change in the number of particles or the radius occurs 
either by capturing one molecule or by emitting another, and since the probability 
of these processes is generally not the same, the fluctuation of n will be subject 
to laws similar to those for the movement of a particle in BROWN’ian motion in 
a force field. The whole problem could be illustrated by the following picture. 
In a vessel, specifically heavier particles are suspended in a liquid. These will 
assume a hypsometric distribution due to the fluctuations in the gravity field*+. A 
certain number of particles is at the height nj — that is the stationary number of 
nuclei. A force is now supposed to act in this cross section that pulls the nuclei 
upwards. The stationary distribution is superimposed by a downward and upward 
diffusion of the particles. The particles that diffuse upwards to large n values give 
the number of droplets that grow into visible structures. Now, in every diffusion 
problem, in order to be able to calculate the number of particles migrating through 
a cross section, the concentration must be given at two points. In our case, the 
downward diffusion must not disturb the stationary distribution, so the equilibrium 
concentration prevails at the bottom, but the concentration of the droplets must 
vanish at the top for very large n. 

We now want to follow what has been said here mathematically and calculate the 
number of droplets with the number of molecules n per cm’. 

The vapor pressure of the droplet is a function of n due to surface tension, and we 


3The idea underlying the strict treatment stems from L. SZILARD. 
4A. EINSTEIN, Ann. d. Physik 19, 371 (1906) 
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will denote it by p(n). Until p(n) is greater than the supersaturation pressure p,, 
the stationary state prevails, i.e. on average as many droplets with the number 
of molecules n will form as will disappear through evaporation during this time. 
This will be the case for all droplets whose particle number is below no. The 
elementary process is always that a droplet captures an impacting molecule or loses 
one through evaporation. 

Let the number of droplets with n molecules be Z(n), then the following equation 
must apply to the stable state under consideration: 


a-p,-Z(n)- F(n) =a-p(nt+1)-Z(n+1)-F(nt+1), 


if a+ p, is the number of molecules, which are separated at 1 cm? per second, a: p(n) 
the number of molecules, which are emitted at 1cm? per second and if ¥(n) or 
F (n+ 1) respectively is the surface of the droplets. We define 
Z'(n) = Z(n)- F(n) 
Z'(nt+1)=Z(n4+1)-A(n+l1). 


For p(n + 1) we can set as a good approximation p(n). We therefore have: 


Z(n+1) _ Pr 


Z(n) p(n)’ 


or, if we make a transformation to In: 


inf2"(n+ I) —InfZ"(n)) _ 
1 =n (5). 


From this equation we obtain a differential equation for Z'(n) if we set 1 = dn. For 
reasonably large n this is justified. The differential equation then reads: 


din(to] ($e) 


dn p(n) 


(1) 


To integrate the differential equation, we need to know p(n) = p(x) as a function 
of n or of the corresponding radius x. 

If we use po. to denote the saturation pressure of a flat surface at temperature T, 
the following relationship applies: 


A&T-\n (=) 220M or Pe) =p; en aT : (2) 
Poo vy 

where M is the molecular weight, y the specific weight of the droplet. This 

derivation assumes that the surface tension o is independent of x, but this is 

certainly not the case if one goes to droplets that are so small that the sphere 

of action of one molecule is not completely filled with others. We will use the 
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equation in this form anyway, which is certainly valid up to about 5 x 1077 cm?, 


but deviations can occur at radii of molecular size. 

p(x) = p, corresponds to a certain radius x = r; this is the radius of the 
characteristic droplet. Its size is determined by the degree of supersaturation. We 
can express the stationary number of nuclei no or generally the particle number n 


by 7 or z: 
4n 5 


i. 
ee cde ee Oy 
gti a (3) 
where L is the LOSCHMIDT number, = 6,06 - 107°, and 
d Ana?y - L 
nm _ Ame*y- Le 4) 
dx M 


Now it is: 


m(g)--m(22)-—m(@)em(B) 


Introducing Equation 2, Equation 3, Equation 4 and Equation 5 in Equation 1 
yields: 


dIn[Z'(n)]__ dIn[Z’(n)}_ MM CD in { 2 
dn _ dx Are?y- Lo Try ee a 


or 


dIn[Z'(n)] 8220 | dre? 2 -@T-In (2) 
de ae kT 
By integration for x we get the function we are looking for: 


4nato + 82-2 BT -In( ) 


De, 
Z'(n) = F(n)-Z(n) =C- exp < — oF ue 


One can also write it a little differently if one introduces the characteristic radius. 


Then 
FSS SURI S665 {a . (1 ss) . 


 3r 
or Z(n) itself, the number of droplet with n molecules per cm?: 


Oe FG ees {3 , (1 = 3 )} | 


The constant in front of the e-power is kept as unknown because one cannot use 
the usual integration because the formula is not valid in a significant area of the 
appropriate range. 


5VAN DER WAALS-KOHNSTAMM, Lehrbuch der Thermodynamik. 
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However, the droplets whose radius have reached r will behave differently than the 
droplets that have a p(n) > p,. By occasionally adsorbing molecules, they will be 
stabilized as their pressure drops below p,, and they will quickly grow into visible 
structures. According to our picture given at the beginning, this corresponds to 
the force that acts on the droplets in the cross section no and pulls them upwards. 
Our Equation 1 will no longer apply to this n region, since there is no longer a 
stationary distribution in the latter. 

a:p,: F(n)- Z(n) will be greater than a- F(n+1)-p(n+1)-Z(n+1) and the 
excess is nothing more than the number of particles migrating through the cross 
section no towards larger n: 


J=Z(n)-prp-a-F(n)—Z(n+1)-a-p(n+1)-F(n+1). 
We set: 
p(n) = p(n +1). 


By J we mean the number of droplets formed per second. We can reformulate this 
equation by adding and subtracting a- p(n): #(n)-Z(n). This yields: 


—J = [Z'(n+1)— 2'(n)] -p(n)-a+ Z"(n)-a- [p(n) — py] . (6) 


We don’t need to take the In now, but the equation already approximately represents 
a differential equation, since Z’(n) does not vary much with n in this range: 


J dZ"(n) 
a-p(n) dn 


If we put for shortness: 


= y(n 
On) cae 
we get the linear differential equation: 
dZ'(n) , Aj 
-Z = 
Fe + elm) 2'(n) + =0, 


which has the solutions: 


a(n) = ex {- [etnyan) [af exp {+ [ atnyan an] | 


y(n) can be approximated for the small range n — np = pw with the following 
functions, namely for positive ju: 


P(r) = Pr _ 2 


y(n) = . —a’-(n— 1) =a? +p, (7) 
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for negative 1: 
y(n) = In (=) (8) 


a- p(n) changes even less with n, since the condensation of a molecule hardly 


changes the surface area and the vapor pressure. So we can take at in front of 


the integral and set p, = p(n). 
In Equation 6 —a? is the slope of the curve at the point n = no or v =r, wp = 0. 
One can calculate —a? from the curve y(n): 


d d (p(n) — py d p(n) 
2 
—a* = | — = 9 eee 
F060) n=no (5, ( Pr n=no dn ‘ Pr n=no 
d 20M d 20M 
ep Prgese ] ~exryrRT _ POF ft etre 
a (= ( N(Deo) “6*7# )) (+ (Sa7)) 
P d 20M dz 
-a = . 
deaty AL din), 
a 20M M 
os ey AT Anny: |. 


We therefore find: 


>  8no0M? 
VF RT-L 
or for a itself after some rearrangement: 
Pr kT 
=ln{—). : 9 
eae ( Be) 20 F (9) 


where .¥ is the surface of the characteristic droplet, —a? is of the order of magnitude 
10~? for water at T’ = 300°C. 
J 


In the solution of the differential equation, after we have taken aT in front of 
a2 2 

the integral, the following integral hs e- 2 dj remains to evaluate. However, we 

can extend this probability integral to oo instead of yz, because the weight of the 


integral is between 0 and wp: 


n a lee) 
/ eS e(njan _ | ge ay? ai | ee a? ‘ine as + STE 
no 0 0 a 2 


and therefore we get: 


and if we put 
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is follows that — 
Z'(n)-e°? =B-J-A. 


To determine the two constants in the solution of the differential equation, we 
introduce the following boundary conditions. For large positive pz Z’(n) must vanish 


a? 2 . 
and even more so Z’(n)-e~ 2, so it follows: 
B-J-A=0. 


On the other hand, for negative 4, Z’(—j) transforms into the equilibrium distri- 
bution. We get, because: 


and therefore: 


But one can proof that: 


namely 


zcw=en{ f"n( Ps) an), o8Paen{ fn (A) a] 


and therefore: 


Thus, we get for J: 
Z' (no) C _ Anr2o 1 C Fo 


J = —— = —_. kT 3 = ——-+€ 3kT 


A. OA 2A , 


or, if we substitute for A its value and since a-p, is the average number of collisions 
per cm? and second according to the kinetic gas theory: 


1 


rr 11 
e e Ir RTM (11) 
we have for J as a final result: 
Pr L F Pr _ Feo 
= l cara 12 
J oF C ys Vea (=) ESRF ol (12) 
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or, after some minor rearrangement: 


Byes 20LM  _ #6 
— Fy&RT T 


J 


-E€ 3kT = 


This is the number of drops formed per second, which remains constant as long as 
the initial conditions are noticeably maintained. It is hereby strictly proven that 
it is proportional to the number of nuclei that would be present in a stationary 
state if no further growth would occur, in accordance with the assumption made by 
VOLMER and WEBER. In addition, it was possible to specify the factor by which 
this number must be multiplied in order to obtain the reaction rate, namely the 


drop formation per second. 
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